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We show that the limiting eigenvalue density of the product of n identically distributed random 
matrices from an isotropic unitary ensemble (IUE) is equal to the eigenvalue density of n-th power of 
a single matrix from this ensemble, in the limit when the size of the matrix tends to infinity. Using 
this observation one can derive the limiting density of the product of n independent identically 
distributed non-hermitian matrices with unitary invariant measures. In this paper we discuss two 
examples: the product of n Girko-Ginibre matrices and the product of n truncated unitary matrices. 
We also provide an evidence that the result holds also for isotropic orthogonal ensembles (IOE). 
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Introduction 

Free probability theory is a fusion of non-commutative probability theory and the concept of free independence. 
Since 1991, when the link between the free probability theory and random matrix theory was established pQ, several 
new results have been proven in an easy and powerful way in the limit of infinitely large random matrices [2H1] ■ In 
this note we demonstrate a simple albeit quite counterintuitive result that the spectral density of the product of n 
free, identically distributed random matrices from an isotropic unitary ensemble (IUE) is equal to the spectral density 
of the rt-th power of a single matrix from this ensemble in the limit of infinite matrix size. The proof is based on the 
multiplicative properties of the S-transform and the Haagerup-Larsen theorem [5]. 

The motivation for the present work comes from the observation made in [5HS] that the eigenvalue density of n 
independent Girko-Ginibre [HI EH] matrices is identical to the eigenvalue density of the n-th power of a single Girko- 
Ginibre matrix in the limit of infinite size. This observation leads to the question whether this is a feature of only this 
particular class of matrices or if there exists a larger class of matrices that have this property. In the present paper we 
show that there is indeed a larger class of matrices sharing this property - a class of random isotropic matrices. We 
begin with defining isotropic matrices. Then we present the main result in detail and its derivation. Finally we outline 
a few sample applications, related to the recent interest in the literature. In particular we apply our result to the 
product of Ginibre-Girko matrices and rederive the density known from J5HE]. We also consider classes of truncated 
unitary and orthogonal matrices and compare our predictions to Monte Carlo simulations, which allow us to identify 
finite size corrections. We conclude the paper with a short summary. 

Isotropic random matrices 

It is convenient to introduce the concept of isotropic random matrices in analogy to isotropic complex random 
variables z that have a circularly symmetric probability distribution depending only on the module \z\. Using polar 
decomposition, one can write z = re 1 ^ where r is a real non-negative random variable and is a random variable 
(phase) with a uniform distribution on [0, 2-k). Isotropic random matrices are defined by a straightforward general- 
ization of isotropic complex random variables. A square N x N matrix x is said to be isotropic random matrix if it 
has a polar decomposition x = hu in which h is a positive semi-definite Hermitian random matrix and it is a unitary 
random matrix independent of h and distributed on the unitary group U{N) with the Haar measure. In short, u is 
a Haar unitary matrix. The random matrix h plays the role of the radial part of x. Such random matrices form an 
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ensemble of isotropic unitary matices (IUE). An example is an ensemble generated by the partition function [TTl I12j : 

Z = J Dx e -^v^ x ) (1) 

where Dx = \~\ i3 d(Rexij)d(lmxij) is a flat measure, and V{a) is a polynomial in a. Another natural class of IUE 
matrices are matrices of the form x = vdu where d is an NxN diagonal matrix having real positive random eigenvalues 
with the given probability distribution and v and u are two independent Haar unitary matrices on the unitary group 
U(N). By analogy one can also consider isotropic orthogonal ensemble (IOE) given by the decomposition x = so with 
s being a positive scmidefinite real symmetric matrix and o being a Haar orthogonal matrix. In this case, when one 
considers an ensemble given by a partition function like that given above, one has to replace a; by a real matrix. 

In mathematical literature, isotropic matrices for N — > oo are called i?-diagonal [13]. In this note we prefer to call 
them isotropic (or IUE, IOE) in the large N limit. IUE matrices have an eigenvalue distribution independent of the 
polar angle on the complex plane. In the limit when the matrix size N — ¥ oo one can find an explicit relation between 
the eigenvalue density of the matrix h 2 and of the matrix x [31 [TTJ [T^l [H] . We briefly recall this relation below. Let 
us mention that the angular independence of the eigenvalue density does not imply that the matrix is isotropic. For 
example a block diagonal matrix of the form: 

where hi, hi are independent Hermitian matrices of dimensions N\ and N2, N1 + N2 — N, and ux, U2 are Haar unitary 
matrices on U(Ni) an U(Nz) respectively, has a circularly symmetric eigenvalue density in the complex plane, but it 
is not isotropic. Intuitively, this is because the split into u\ and U2 breaks the isotropy in the whole U(N) group. 

Main result 

The main result of this paper is as follows: consider n identically distributed isotropic matrices x\, x%, . . ., x n 
generated independently from a given IUE (isotropic unitary ensemble). In the limit N — > 00 the eigenvalue density 
of the product X n = x\Xi ... x n becomes identical as the eigenvalue density of the n-th power x n of a single matrix 
x from this ensemble (e.g. x = x{). In other words, the probability that a randomly chosen eigenvalue of X n lies 
within a circle of radius r: Prob(Ax„ < r) approaches for A — > 00 the probability that a randomly chosen eigenvalue 
of x" lies within the same circle: Prob(A" < r). One can use this observation to derive the eigenvalue density of the 
product X n = X1X2 ■ ■ ■ x n if the eigenvalue density of x is known. In particular one can immediately show that the 
eigenvalue distribution of the product of n independent Girko-Ginibre matrices has a simple form: 

p(z,z) = — \z\- 2+2 ' n for \z\ < 1 (3) 
7rn 

and zero for \z\ > 1, in agreement with ITS1 IT^l l2"2l |2"3"] . It is interesting to note that the matrices X n X^ obtained 
from the products X n of Girko-Ginibre matrices generate a Fuss-Catalan family of distributions [15] that have however 
a much more complicated limiting eigenvalue density [16) . Another interesting case is the product of n independent 
truncated unitary matrices [T7j that is 

p(z,z) = ^|zr 2 + 2 /"(l-M 2A r 2 for \z\<(-^-) n ' 2 (4) 

717T \ 1 + K ) 

and zero otherwise. The truncated matrices have dimensions N x N. They are obtained by removing L columns and 
L rows from the original (N + L) x (N + L) Haar unitary random matrix. The result holds for N — > 00 and k = L/N 
fixed. 

This is a counterintuitive result, so let us stress that it only holds in the limit N — > 00. For finite N the eigenvalue 
distributions of the product of X\ . . . x n and of the power x n differ. The difference however disappears when N tends 
to infinity, as we illustrate it below. 

Derivation 

Consider an IUE ensemble of random matrices x = hu of dimensions N x N . In the large N limit the random 
matrices can be represented as free random variables and one can use the Haagerup-Larsen theorem [S] that relates 
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the eigenvalue density of x to the eigenvalue density of h 2 by the following formula: 

S h , (F x (r) - 1) = i . (5) 

where F x (r) is the cumulative density function for the density of eigenvalues of x on the complex plane and Sy^ (x) is 
the S-transform for the matrix h 2 . The cumulative density function 

F x (r)= / d 2 zp x (z,z) = 2tt / dssg x (s) = / dsp x (s), (6) 

J|z|<r JO JO 

can be interpreted as the fraction of eigenvalues of x in the circle of radius r centered at the origin of the complex 
plane. It is related to the eigenvalue density p x (z,z) = g x {\z\) that depends on the distance from the origin \z\. The 
integrand dsp x {s) — 2ndssg x (s) is interpreted as the probability of finding eigenvalues of x in a narrow ring of radii 
\z\ and \z\ + d\z\: 

F' x {r) = Px (r) = 2nrg x (r) . (7) 

The prime denotes the derivation with respect to the radial variable. The cumulative density function F x (r) enters 
equation ^ as an argument of the S-transform (z) that is related to the eigenvalue density py^ (A) of the matrix 
h 2 (see Appendix A) . The Haagerup-Larsen theorem states also [Til Ell Q2| that the support of the eigenvalue 
density of a; is a ring of radii i? m i n and i? m ax or a disk (if R mm = 0): 

i?Lx= / d\\p h 2(\) , R m 2 n = dAA-W(A). (8) 
Jo Jo 

Let us make few comments. For an i?-diagonal (isotropic) matrix x given by the radial decomposition x — hu, where 
h is Hcrmitian and u is a Haar unitary matrix, the two matrices xx^ = h 2 and x^x = u*h 2 u have identical eigenvalues 
and therefore the S-transforms for xx^ and x*x are identical: S xx f(z) — S x i x (z) = Sf- t 2(z). This means that ^ can 
be written as 

S x t x (F x (r) - 1) = i . (9) 

Let us now apply this equation to the product of n identically distributed i?-diagonal (isotropic) matrices X n = 
x\...x n . The resulting matrix has an identical eigenvalues as H n u n , where H 2 = X^X n so we can apply ^ 
replacing in this equation x by X n : 

SxtxjFxn(r)-l) = ±. (10) 

The S-transform for the matrix X^X n which appears in the last equation can be substituted by the S-transforms for 
individual terms in the product. Indeed, writing 

X\X n = xlxl^X^Xn (11) 

where A„_i = x\ . . . x n -\ we see that 

S xtx n = S Xl_ 1 X n - 1 S xtx n ( 12 ) 

since due to the cyclic properties of trace the moments of xJ l A" ? ' l _ 1 X„_ 1 a;„ are identical as those of x n xl l Xl l _ 1 X n _i 
and moments of XnX^ as those of x^Xn. Applying the last equation recursively we eventually obtain 



S xlx n =I[S xlxt . (13) 



i=\ 



Taking into account that all Xi are identically distributed and having the same S-transform (that we denote by S x i x ) 
we can write the last equation as 

S xlx n = S xix ■ ( 14 ) 
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Inserting this into ( 10 ) we have 



S x , x {F Xn {r)-l) = ^ Jn . (15) 

This equation has an identical form as Q except that on the left hand side F x {r) is replaced by Fx n (r) and on the 
right hand side r is replaced by r 1 /". From this observation it immediately follows that 

Fx n {r) = F x {r 1 l n ) = F xn (r) . (16) 

The last equality follows from the fact that eigenvalues of the matrix x n are equal to the n-th power of the correspond- 
ing eigenvalues of x: F x ™{r) = Prob(|A|™ < r) — Prob(|A| < r 1 /™) = F^r 1 /"). So we see that indeed the product of n 
identically distributed isotropic matrices X n — X\x 2 ■ ■ ■ x n has the same eigenvalue distribution as the n-th power x n 
of a single matrix in the product. In practice, the eigenvalue distribution of X n can be calculated directly from the 
eigenvalue distribution of a single matrix x by substituting r — > r 1 /™ in the cumulative distribution function F x (r) 
(JGh . The corresponding eigenvalue densities may be found using Q. They read 

px n {r) = -r 1/n - 1 p x {r 1/n ) (17) 



and 



Qx n (r) = -r 2/n - 2 Q x (r 1/n ) ■ (18) 
n 



Applications 

Let us apply these formulas to a couple of examples. First consider Girko-Ginibre matrices [9j [10] that have a 
uniform distribution g x (r) — 1/ir inside the unit circle \z\ < 1. We have 

F x (r) = 2 [ rdr=r 2 for r < 1 (19) 



and 1 otherwise. For the product of n-independent Girko-Ginibre matrices we have ( 16 ) 

F Xn (r) = r 2/n for r < 1 (20) 
and one otherwise. Taking the derivative with respect to r Q we find the corresponding densities: 

Px n (r) = -r 2/n - 1 6(l-r) (21) 
n 

and 

QxAr) = — r 2 /"- 2 #(l-r) (22) 

where denotes the Heaviside step function. This result agrees with that obtained using different methods in 
[51 HH1 HH1 1221 123] as mentioned in the introduction of the paper. 

As the second example, we consider the product of n truncated unitary matrices |17j . The cumulative eigenvalue 
distribution of a single matrix from this ensemble is 

w = t^ fOT ^(rb) V2 ( 23 ) 

and 1 otherwise. The coefficient n — L/N is the ratio of the number of rows and columns L removed from a Haar 
unitary matrix of dimensions (N + L) x (N + L). This truncation leaves a matrix of dimensions N x N. In Appendix 
B we show how to derive this result using free random variables. The corresponding density reads: 

k ( ( 1 x 



to (r) = -(l-rT'*n — ) -r) . (24) 
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FIG. 1: (Color online) Numerical verification of theoretical formulas (J3J (a) and B (b),(c) for the radial part p(r) = F'(r) 
of the mean spectral density p(z,z) of the product of independent matrices, (a) Numerical histograms for the product of 
3 independent Gaussian random matrices TV = 400 (black crosses), TV = 200 (red circles) and TV = 100 (blue rotated crosses) 
compared to theoretical prediction for TV — > oo (solid green line). Each histogram is made for 10 7 eigenvalues. The numerical 
histograms approach theoretical curve as the size of matrices increases, (b) An analogous plot to (a) for the product of 2 
independent truncated unitary matrices with ratio k = | and TV = 360 (black crosses), TV = 180 (red circles) and TV = 90 (blue 
rotated crosses). Each histogram is made for 9 x 10 6 eigenvalues, (c) An analogous plot to (a) and (b) for the product of 3 
independent truncated unitary matrices with ratio k = \ and TV = 320 (black crosses), TV = 160 (red circles) and TV = 80 (blue 
rotated crosses). Each histogram represents 8 x 10 6 eigenvalues. 



Using ( 16 1 we find the distribution of eigenvalues for the product of n such matrices: 



r 2/n / l \ n/2 



F ^{r) = K Y — 2 j^ for r<l — ) (25) 



and 1 otherwise. The corresponding eigenvalue density is 



9x^r)^^r^-\l-r^)-^((-^-X' 2 - r) . (26) 

7T-7T \ \ 1 + K J I 



Numerical comparison and finite size effects 



In order to crosscheck our results, we use Monte-Carlo simulations for generating (sampling) finite size random 
matrices from ensembles in question. An agreement between the analytical formulas ^ or Q and numerical results 
is observed taking into account finite size corrections. The shape of obtained distributions p(r) — F'(r) ([7| is shown in 
the figure [TJ In the TV — > oo limit, distributions have got compact support, and the sharp drop at the edge is present. 
For finite TV the spectra do not have a sharp threshold - instead they tend to zero continuously in an extended 
crossover region, and the difference between the product of independent matrices and the corresponding power of a 
single one is visible in this region (figure [2]). The eigenvalue density of the product of independent matrices approaches 
the theoretical curve faster than of the corresponding power of a single matrix. Only radial distributions p(r) — F'(r) 
([7]) are shown, since eigenvalue densities are circularly symmetric on the complex plane. The shape of the finite size 
corrections for Girko-Ginibre distribution was discussed in [TBI 119] , 

We also performed numerical simulations for the products of truncated orthogonal matrices as an example of 
multiplication of IOE matrices. In the large TV limit both the IUE and IOE densities are expected to have the same 
limiting distribution while for finite TV the distribution in the IOE case is expected to display a characteristic pattern 
that weakly breaks the circular symmetry of the eigenvalue distribution on the complex plane. More precisely, one 
expects that a fraction of eigenvalues accumulates on the real axis and disappears from a narrow depletion region 
close to the axis. The effect was first discussed for real Girko-Ginibre matrices |5D] and later also for orthogonal 
truncated matrices [21 . It is known to be a finite size effect in the sense that the fraction of eigenvalues forming the 
pattern tends to zero for TV — > 00 so the full circular symmetry of the eigenvalue density is restored in the limit. In 
fact, this is exactly what we see in our numerical simulations of the product of truncated matrices. First we observe 
that the radial distribution of eigenvalues of product of two truncated unitary matrices is identical to the case of 
truncated orthogonal matrices except in a small region close to r = (see fig. [3ja). In figure [3jb we compare finite 
size distributions for the product of IUE (lower part) and IOE (upper part). We see that the IUE distribution is 
circularly symmetric up to the statistical noise while the IOE distribution has an elongated shape close to the real 
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FIG. 2: (Color online) Numerical analysis of the finite size effects for the radial part p(r) = F'(r) (TTb of the mean spectral 
density p (z, z) of the product of independent matrices in comparison to the power of a single matrix, (a) Numerical histograms 
for product of 3 independent Gaussian random matrices TV = 200 (black crosses) and one matrix raised to 3'rd power for 
TV — 200 (blue circles) compared to theoretical prediction for TV — > oo (solid green line). Each histogram is made for 10 7 
eigenvalues. Plots are zoomed in the region, where the difference in the shape is visible, (b) An analogous plot to (a) for the 
product of 2 independent truncated unitary matrices (black crosses) and 2'nd power of a single truncated unitary matrix (blue 
circles) with ratio k = | and TV = 180. Each histogram is made for 9 x 10 6 eigenvalues, (c) An analogous plot to (a) and 
(b) but the product of 3 independent truncated unitary matrices (black crosses) and 3'rd power of a single truncated unitary 
matrix (blue circles) with ratio k = \ and TV = 160. Each histogram is made for 8 x 10 6 eigenvalues. 




FIG. 3: (Color online) Numerical comparison of the eigenvalues of product of two truncated orthogonal and unitary matrices, 
(a) The radial part p(r) = F'{r) |7| of the mean spectral density p(z,z) for unitary (black crosses) and orthogonal (red circles) 
matrices for TV = 100, k = |. Each histogram is made for 9 x 10 6 eigenvalues. The theoretical prediction for TV — > oo is shown 
for comparison (solid green line), (b) Full eigenvalue distribution of orthogonal (upper half of complex plane) and unitary 
(lower half of complex plane) truncated matrices for TV = 100, k = |. (c) Full eigenvalue distribution of orthogonal truncated 
matrices for same TV,k parameters. The real eigenvalue band is clearly visible. 

axis, as expected. Finally in fig. [3jc we show the full spectrum on which one can clearly see an accumulation of 
eigenvalues on the real axis. 

Discussion 

In this note we have shown a simple, (and as far as we know) new relation between the spectral properties of the 
product of n identically distributed isotropic random matrices from the given IUE ensemble and spectral properties 
of n-th power of a single matrix from this ensemble. We stress a nonintuitive aspect of this result that tells us 
that independent matrices, when multiplied, give the same eigenvalue density as the product of fully correlated 
(identical) matrices. In a sense it is a self-averaging effect: a single random matrix from a isotropic ensemble is a 
good representative to describe products of matrices from this ensemble in the limit TV — ¥ oo. - 

We have supplemented our analytic proof with extensive numerical simulations, allowing us to see how the finite 
size effects vanish in the thermodynamical limit. For Girko-Ginibre finite size effects agree with those conjectured in 

Ha nsi. 
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Our result elucidates the transparent analytic structure noted in several recently published papers on the products 
of random matrices [6HS1 HTKT91 [221429] and provides a powerful tool for the derivation of similar results for products 
of some application-designed, isotropic random matrices of large (infinite) size. 
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Appendix A 



In this appendix we briefly recall basic facts about the S-transform, introduced by Voiculescu in free random 
probability |30j . Consider a Hermitian random matrix a. One usually defines the Green function 

G a (z) = lim 1 (Tr(z - a)" 1 ) = / dX^- . (27) 

N— >oo iV J Z — A 

that is directly related to the eigenvalue density p a (A) . Note that the density is a function of a real variable while the 
Green's function is a function of a complex variable. The Green's function generates moments pb a k (if they exist) 

G ^)='+E|^ (28) 
fe=i 

of the eigenvalue density 

Mafc = lim i (Tr a k ) = f d\p a (X)X k . (29) 

Sometimes it is more convenient to use another generating function, given by a power series in z rather than in 1/z: 

i /i \ 00 
M*) = -Ga( "J -1 = (30) 
^ Z ' fc=l 

and to introduce its functional inverse % a : 

Xa(<f>a(z)) = 4>a{Xa{z)) = Z (31) 

which can also be expressed as a power series in z if the first moment is nonzero: /u i 7^ 0. The S-transform for the 
matrix a is related to the ^-transform as 

Sa(z) = ^Xa(^) • (32) 

The relevance of the S-transform in free probability is related to the fact that it allows one to concisely formulate 
the law of free multiplication. The S-transform of the product of two free (independent) matrices from invariant 
ensembles is a product of the S-transforms of individual matrices: 

Sab(z) = S a (z)S b (z) (33) 

The multiplication law was formulated in free random probability 30J but it also can be rederived in random matrix 
set-up using field theoretical techniques for the summation of planar Feynman diagrams and it can be generalized to 
non- hermitian matrices |2"2"1. 



Appendix B 



In this appendix we rederive the distribution of a single unitary truncated matrix ( 23 ) using free probability and 
the Haagerup-Larsen theorem. We first construct the density of an (N + L) x (N + Ljmatrix y = pu where 

p = diag(l, . . . , 1, 1 __0) (34) 

N L 

is a projection matrix and u is an Haar unitary matrix of dimensions (N + L) x (N + L). In order to calculate 
the S-transform fo r th e projector p = p 2 we first observe th at a ll moments of p are equal ^ = N/(N + L). Hence 
= ivfri=5 < 301 ' *p = n/(nIl)+z < 31 > and eventually (32 1 



1 + z 



N/(N + L) + z 



Inserting this to ^ we find 



F y (r) 



L 

n + lY 



for r < ^N/{N + L) 



(35) 



(36) 



and 1 otherwise. We see that F y (Q) — L/(L + N). This means that there are L eigenvalues equal zero. They are 
inherited from the zero eigenvalues of the projector. We can now reduce dimensionality of the matrix y by removing L 
zero eigenvectors. The remaining matrix x that has no trivial zero eigenvalues. This gives the result given in equation 
( 23 1 for a single truncated matrix 



F x {r) 



N + L 
N 



FJr) 



L 



N + L 



Nl-r 2 



for r < y/N/(N + L) 



(37) 



and one otherwise. The term —L/(N + L) removes L zero eigenvalues out of (N + L) eigenvalues, and the factor 
[N + L)/N restores the total normalization. 
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